Polar behavior of double perovskite (Bi,Pb)ZnNbOe and (Bi,Sr)ZnNbOe: Density 

functional calculations 
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The polar behavior of double perovskite (Bi,Pb)ZnNbOe and (Bi,Sr)ZnNb06 was investigated 
using first principles density functional calculations within the local density approximation. These 
materials have both A-site size disorder and ions with sterochemical activity. We found a strong 
ferroelectric distortion in both materials. The polarization of the Pb material is >80 (iC/cm 2 
along the rhombohedral direction and the Sr based material is only slightly inferior. The high 
polarization results mainly from a large off-centering of Bi and the large Born effective charge of 
Nb, which yields a large contribution polarization although the Nb off-centers by a smaller amount. 
Neither of these materials favors a tetragonal state, and therefore solid solutions with PbTiOs may 
show morphotoropic phase boundaries. 

PACS numbers: 77.84.Ek,71.20.Ps 
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I. INTRODUCTION 

Perovskite ferroelectrics and relaxor ferroelectrics are 
widely used in applications including electromechanical 
transducers, actuators, passive electronic components, 
memory devices and microwave technology. There is 
strong interest in obtaining improvements in materials 
performance, such as high polarization, high Curie tem- 
perature, improved electromechanical coupling and lower 
loss. One approach is through the exploration of new ma- 
terials. Besides this technological motivation, these ma- 
terials have been the focus of considerable fundamental 
interest, particularly in unraveling the complex interplay 
between crystal chemistry and lattice distortions and in 
understanding the range of polar behaviors that can be 
produced by strain and chemical substitutions. Much has 
been learned over the past several years. 

The lattice instabilities of perovskite oxides, ABO3 
are often understood using the tolerance factor t — 
(ro+ryi)/ V2{ro+ r B), where ro, ta and tb are the ionic 
radii for the O, ^4-site and B-site ions, respectively^— 
These materials are classed in two groups according to 
whether the tolerance factor is less than unity or not. 
Several t > 1 ferroelectrics are known, including the pro- 
totypical materials, BaTiC>3 and KNbC>3. In these mate- 
rials there is a tendency towards lattice distortions where 
the too-small B-site ion off-centers. This is strongly en- 
hanced by covalency involving O 2p states and nomi- 
nally unoccupied -B-site d states. This then leads to 
enhanced Born effective charges and large longitudinal 
optic - transverse optic splittings. 4 ' 5 In the majority of 
t < 1 materials the BOq octahedra are tilted and they 
are not ferroelectrics^ The exceptions are generally per- 
ovskites with lone-pair A-sites such as Pb and Bi. In 
these cases, the ideal cubic perovskite structure is often 
unstable against octahedral tilting, but the ferroelectric 
instability, which is enhanced by the stereochemical ac- 
tivity of the A-site, is stronger. In materials with t < 1 
but no stereochemically active A-site ion, there is often 



a substantial ferroelectric instability of the ideal cubic 
structure, but now the tilt instability is stronger, and 
so ferroelectricity does not occur X"— However, accord- 
ing to first principles calculations, if the octahedra are 
prevented from tilting by mixing large and small A-site 
ions they may become ferroelectrics. 3 There is exper- 
imental support for this. For example, ferroelectricity 
in the (Ba,Ca)Ti03 solid solution is strongly enhanced 
compared to what would be expected interpolating be- 
tween the properties of the end-point compoundsfi 1 - Im- 
portantly, this leads to new compositions that have high 
piezoelectric coefficients i 11 ! 12 

Generally, high piezoelectric coefficients and elec- 
tromechanical coupling are found near phase boundaries. 
In particular, the electromechanical coupling in per- 
ovskite piezoelectrics is associated with polarization ro- 
tation near morphotropic phase boundaries (MPBs)j 13 ' 14 
These are generally boundaries between ferroelectric 
states differing in the direction of their ground state 
polarization. Therefore one strategy for finding new 
piezoelectric compositions is to search for solid solutions 
with MPBs. A particular emphasis has been in solid 
solutions with PbTi03 (PT) as one end-point, follow- 
ing the commonly used Pb(Zr,Ti)03 (PZT) system and 
the PbMg 1/3 Nb 2/ 3C>3 (PMN) - PT and PbZn 1/3 Nb 2/3 03 
(PZN) - PT piezocrystal systems^ PT is tetragonal with 
a substantial tetragonality. There is considerable inter- 
est in finding materials that can be alloyed with PT 
and which have strong ferroelectricity with rhombohe- 
dral or related ground states to produce an MPB perhaps 
with high Curie temperature, high polarization and high 
tetragonality on the tetragonal side of the MPB. 

A development of potential importance is the realiza- 
tion that Bi based perovskites can be very useful ferro- 
electrics, as exemplified by progress in BiFeOs. In partic- 
ular, the small ionic radius and stereochemical activity 
of Bi 3+ , leads to large displacements in a low volume 
unit cell, and its high charge (3+ vs. 2+ for Pb) and 
hybridization of Bi 6p states with O 2p states leads to a 
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high Born charged Furthermore, experimental work on 
the BiScC>3 - PT solid solution shows an MPB with en- 
hanced Curie temperature,— and there have been recent 
studies of several other Bi-based perovskites in solid so- 
lutions with PTiiSr— . One problem has been that while 
the small Bi-ion tends to enhance tetragonality on the 
tetragonal side of the MPB, this same effect may prevent 
the occurrence of an MPB, leading to phase diagrams 
with only tetragonal phases. 

Recently, an extremely large tetragonality (c/a ~ 
1.11 and high Curie temperature Tq of ~ 700 °C at 
x = 0.6 was reported for perovskite solution (1 — 
x)Bi(Zn 1/2 Ti 1/2 )03-(x)PbTi0 3 ^ This is reminiscent of 
the behavior of CdTi03 - PT alloys, which also have 
highly enhanced tetragonalityi^ir— In perovskites the 
balance between tetragonal and rhombohedral states is 
generally controlled by a balance between the energy low- 
ering due to tetragonal strain, which favors the tetragonal 
state, and the i?-site off-centering, which favors a rhom- 
bohedral stated (note that with small A-site ions there is 
an additional stabilization of the tetragonal state because 
off-centerings along [001] directions are towards the most 
open face of the O cage around the A-site; therefore, this 
direction is favored for the A-site and this can be impor- 
tant if the displacement is very large^). The interplay 
of A-site and _B-site off-centering, and the importance of 
the _B-site in determining the direction of the ferroelectric 
polarization suggests the exploration of alternate i?-site 
ions in solid solutions with PT to obtain an MPBs. 

Here, we present a study of the polar behavior of per- 
ovskite (Bi,Pb)(Zn,Nb)0 6 and (Bi,Sr)(Zn,Nb)0 6 using 
density functional supercell calculations. The purpose is 
to explore the effect of combining the mechanisms dis- 
cussed above in a perovskite that may be amenable to 
experimental synthesis. Specifically, the compositions 
explored have different size A-site cations (Pb and Bi), 
stereochemically active ions on the A-site (particularly, 
Bi, though Pb may also contribute as discussed below), 
and ions with electronic structures that favor cross-gap 
hybridization and therefore ferroelectricity on the i?-sitc 
(Nb and Zn). 

As mentioned, one motivation for this study is to ex- 
plore the use of different size A-site ions in relation to 
stereochemical activity to obtain A-site driven ferroelec- 
tricity. The ionic radii 2 - of the A-site ions considered here 
are r Pb 2+ = 1.63 A, r Sr 2+ = 1.58 A, and r Bi 3+ = 1.31 A. 
Thus, Bi 3+ is significantly smaller than the other A-site 
ions and Pb and Sr have approximately the same size, 
but unlike Sr, Pb has stereochemical activity that can 
favor ferroelectricity, as in PbTi03. 

Related to this there is an experimental re- 
port of the synthesis and some physical prop- 
erties of perovskite (BiQ.5Sro.5)(MgQ.5Nbo.5)03, 
(Bio.5Bao.5)(Mgo.5Nb .5)0 3 , (Bio.5Sro.5)(Zno. 5 Nbo.5)0 3 , 
(Bi . 5Bao.5)(Zn .5Nbo. 5)03, and the solid solu- 
tion of (Bi sBao 5)(Zn sNbo 5)03 with PT and 
(Bio.5Sro.5)(Mgo.5Nbo.5)0 3 with SrTi0 3 ^^ The com- 
pound (Bi . 5Sr .5)(Mg .5Nbo. 5)03 was reported to be 



ferroelectric based on the observation of a hysteresis loop 
below 103 K. All four compounds were reported to show 
phase transitions above room temperature and relatively 
high dielectric constants. They also show substantial 
microwave loss, which may indicate ferroelectricity. 



II. APPROACH 

The perovskite alloys were studied using supercell 
calculations within density functional theory. The 
main results are for compositions BiPbNbZnOe and 
BiSrNbZnOg. The large charge difference between Zn 2+ 
and Nb 5+ may be expected to lead to a strong ordering 
tendency on the B-site. This is the case in the relaxor 
systems PbZn 1/3 Nb 2 /303 (PZN) and PbMg 1/3 Nb 2 / 3 03 
(PMN), where a double perovskite like structure occurs 
even though the Nb:Zn stoichiometry is not 1:1 (specif- 
ically, in PZN and PMN one double perovskite B-site 
sublattice is Nb and the other is 1/3 Nb and 2/3 Zn or 
Mg, which in turn may further order) j^2t— In the present 
case the stoichiometry is 1:1, which should further stabi- 
lize the ordering. Therefore we assume an ordering of the 
-B-site lattice into a double perovskite structure. In con- 
trast, there is no a priori reason to assume that the A-site 
lattice consisting of Pb and Bi is chemically ordered. 

BiPbZnNb0 6 (BPZN) and BiSrZnNb0 6 (BSZN) have 
ferroelectrically active ions on both the A and B sites. 
On the B-site both Zn 2+ and Nb 5+ favor ferroelectricity, 
as is evident from the stronger relaxor ferroelectricity of 
PZN relative to PMN and the ferroelectricity of KNb0 3 . 
Likewise Bi 3+ and Pb 2+ are stereochemically active, and 
favor ferroelectricity on the perovskite A-site. However, 
both BPZN and BSZN are expected to be strongly A- 
site driven materials both based on ionic radius (t < 1, 
and A-site size disorder) and lone pair physics ^2 and as 
such a competition between octahedral tilts and polar off- 
centering may be expected similar to the Pb(Zr,Ti)03 
system. 33 Therefore it is important that the supercells 
are constructed in a way that this competition can be in- 
cluded. The minimum size for this is 2x2x2, or 40 atoms, 
since such a cell has an even number of units along the 
[001], [011], and [111] directions and therefore can accom- 
modate the various Glazer tilt patterns 

An additional constraint comes from the fact that with 
multiple A-site ions it is possible to choose cells that have 
polar spacegroups just because of the cation ordering and 
not because of lattice instability. This would not repre- 
sent ferroelectricity, but would rather simply be an arti- 
fact of the selected order. Therefore even though Bi and 
Pb are not expected to chemically order in the solid solu- 
tions, we consider them in a highly ordered state within 
our supercells. Here we select a rock-salt ordering of 
Pb and Bi on the A-site. The cation ordering with this 
choice has symmetry F-i3m, which is non-polar. Thus 
any polarization or off-centering is a consequence of lat- 
tice instability and not the choice of cation ordering, and 
furthermore, one may expect that a more disordered alloy 
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would have a better frustration of the tilts and a stronger 
ferroelectric tendency^ 

The density functional (DFT) calculations were done 
using two methods. We used the general potential lin- 
earized augmented planewave (LAPW) method^ for the 
structure relaxations. This is a full potential all electron 
method. We did the calculations using the local density 
approximation at a scalar relativistic level. Well con- 
verged basis sets were employed with LAPW sphere radii 
of 2.3 Bohr for Pb, Sr and Bi, 1.90 Bohr for Zn and Nb 
and 1.55 Bohr for O. We included local orbitals to relax 
linearization and to accurately treat semicore states^ 
The polarization and Born effective charges were calcu- 
lated with an ultrasoft pseudopotential method as im- 
plemented in the quantum espresso package^ using the 
relaxed structures from the LAPW calculations. The re- 
sults were converged with respect to the Brillouin zone 
sampling, which was tested. The polarization was ob- 
tained using the Berry's phase method. 



III. RESULTS 

We started by assuming a pseudocubic structure, i.e. 
setting the lattice parameters of our 40 atom super- 
cells to be orthogonal and equal. We then fully relaxed 
the internal coordinates of all atoms in the cell, with 
no symmetry constraints. This was done as a function 
of the lattice parameter to obtain the equilibrium cell 
volumes for BiPbZnNb0 6 and BiSrZnNb0 6 . This was 
done in the LDA using the LAPW method. The ef- 
fective perovskite lattice parameters for the minimum 
energy were 3.99 A, for BiPbZnNb0 6 and 3.97 A, for 
BiSrZnNb0 6 . As mentioned, (Bi .5Sro.5)(Zn .5Nbo.5)0 3 
was synthesized by Kosyachenko and co-workers j2£ They 
reported a weakly tetragonal structure with lattice pa- 
rameters a=&=4.002, and c=4.015 A (c/o=1.003) at 293 
K. These values are ~ 0.9 % larger than our calculated 
LDA lattice parameter. This size of underestimation is 
typical of LDA errors. 

The structure of the 40 atom pseudocubic 
BiPbZnNbC>6 supercell at the LDA equilibrium volume 
is depicted in Fig. [TJ As may be seen, there are large 
off-centerings of the Bi ions and also some octahedral 
tilting. Fig. [5] shows the cation off-centerings in this cell 
for the different sites with respect to the centers of their 
O cages (the 12 nearest O ions for the A-site and the 
6 nearest for the £?-site). The largest off-centerings are 
of the Bi ions, consistent with Fig. [TJ and what might 
be expected based on the small size and stereochemical 
activity of Bi 3+ . These displacements are large (0.64 
A) and collinear along a [111] direction. Importantly, 
there are also large collinear displacements of the other 
cations. The Pb ions displace by ~ 0.39 A, the Zn by ~ 
0.27 A, and the Nb by ~ 0.27 A, all very close to [111]. 
This cooperative nature of the ferroelectricity, where 
all cations off-center significantly, is a characteristic of 
good perovskite ferroelectrics such as KNbC>3, BaTiC>3 
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FIG. 1: (color online) Structure of the 40 atom pseudocu- 
bic BiPbZnNbOe supercell after relaxation. The B-site - O 
bonds are shown. The O ions are shown by (dark) red, Zn by 
(medium) gray, Nb by smaller (medium) blue, Pb by gray and 
Bi by (smaller, light) pink spheres. For clarity eight (2x2x2) 
supercells are shown. The top panel shows a view along (001) 
and the bottom shows a view along (011). Note that in addi- 
tion to the cation off-centerings, there are noticeable octahe- 
dral tilts. 




Bi Bi Bi Bi Pb Pb Pb Pb Zn Zn Zn Zn Nb Nb Nb Nb 

FIG. 2: (color online) Cation off-centerings along the Carte- 
sian directions with respect to their O cages for the various 
sites in LDA relaxed structure of the pseudocubic 40 atom 
supercell of BiPbZnNbOe at a lattice parameter of 3.99 A. 
The total displacement magnitude of each ion is also given. 
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FIG. 3: (color online) Electronic DOS of the pseudocubic re- 
laxed BiPbZnNbOe supercell on a per formula unit basis. The 
top panel shows the total DOS and the Bi s and Pb s projec- 
tions. The bottom panel shows Bi p, Pb p, Nb d, Zn sp and 
O p projections. The projections are the integrals of com- 
ponents of the charge density with given angular character 
within the LAPW spheres. Note that for extended orbitals 
this gives an underestimate. The energy zero is set at the 
valence band maximum. 



and PbTi03 i 27 ' 38 These displacements would yield a 
point charge polarization based on nominal charges 
of 58 ^(C/cm 2 using nominal charges, with the largest 
contributions coming from the Bi and Nb (note the +5 
charge state of Nb). The actual polarization is higher 
(see below) because the nominal charges are enhanced 
due to covalency. 

As is usual in perovskite ferroelectrics, we find that the 
ferroelectricity is enhanced as the volume is expanded 
and conversely under compression. However, over the 
volume range investigated (pseudocubic lattice parame- 
ter, 3.88 A < a < 4.07 A), the material remains fer- 
roelectric, with cation displacements along [111]. The 
average cation displacements as a function of a are given 
in Table HI The LDA typically underestimates lattice pa- 
rameters, often by ~ 1-2 %, and as such the prediction of 
ferroelectricity in this material would seem to be robust 
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FIG. 4: (color online) Cation off-centerings along the Carte- 
sian directions with respect to their O cages for the vari- 
ous sites in LDA relaxed structure of the 40 atom supercell 
of BiPbZnNbOe at at the same volume as in Fig. [21 but 
with an imposed c/a ratio of 1.0606. Note that despite the 
large tetragonal strain the cation displacements remain nearly 
collinear and do remain reasonably close to a [111] direction. 



TABLE I: Average displacements of the cations relative to 
their O cages in relaxed 40 atom BiPbZnNbOe cells as a func- 
tion of pseudocubic lattice parameter, a. 



a(A) 


MA) 


<MA) 


<MA) 


^Nb(A) 


3.88 


0.50 


0.28 


0.20 


0.22 


3.92 


0.54 


0.31 


0.23 


0.24 


3.98 


0.62 


0.38 


0.26 


0.27 


3.99 


0.64 


0.39 


0.27 


0.27 


4.03 


0.69 


0.44 


0.30 


0.30 


4.07 


0.75 


0.50 


0.35 


0.32 



against LDA volume errors. 

We now discuss the electronic structure as it relates to 
ferroelectricity. The electronic density of states (DOS) 
and projections for the relaxed supercell are shown in 
Fig. [3 The calculated LDA band gap is 2.7 eV. This is 
almost certainly an underestimate due to the LDA band 
gap error. Considering the sizable gap and the fact that 
all the ions are in chemically stable valence states it is 
likely that BiPbZnNbOe can be made as a good insulator 
at ambient temperature. This is an important consider- 
ation for a ferroelectric material to be used in applica- 
tions. The stereochemical activity of Bi 3+ and Pb 2+ is 
associated with the lone pair chemistry, which is usually 
discussed in terms of a high polarizability of the elec- 
trons in the occupied 6s orbitals of these ions. As shown, 
the Bi 6s states give rise to the DOS peak centered at 
-9.5 eV relative to the valence band maximum (VBM), 
while the Pb 6s states are at -7.4 eV. The O 2p bands 
provide the DOS from -6 eV to the VBM, while the con- 
duction bands are derived from unoccupied orbitals of 
metal character, primarily Nb 4<i, Bi 6p, Pb 6p and at 
higher energies, Zn sp character. There is a modest hy- 
bridization between the Bi and Pb 6s states and the O 
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FIG. 5: (color online) Cation off-centerings along the Carte- 
sian directions with respect to their O cages for the various 
sites in LDA relaxed structure of the pseudocubic 40 atom 
supercell of BiSrZnNbOe at a lattice parameter of 3.97 A. 
The total displacement magnitude of each ion is also given. 
Note the reduced off-centering of Sr relative to Pb (Fig. [2]) 
in BiPbZnNbOe, while the off-centerings of the other cations 
are similar. 



2p states as may be seen from the Bi and Pb s character 
at the top of the O 2p bands. However, this is a mixing 
of occupied states and as such does not signify chemi- 
cal bonding. Instead, as has been noted previously for 
ferroelectric perovskites*^ the important hybridization 
is between the O 2p states and unoccupied metal states, 
which favors ferroelectricity due to chemical bonding and 
leads to enhanced Born effective charges. The energy 
position of the Nb Ad states, which are dominant in the 
lower conduction bands from the conduction band mini- 
mum (CBM) to ~ 1.2 eV above it is important because 
the low energy favors hybridization with O 2p states and 
ferroelectric off-centering, as was discussed in the case 
of PMN and KNb0 3 in relation to KTa0 3 ^^ In the 
present case, there is a strong mixing of O 2p and Nb 
4d states. This is evident from the Nb d character in 
the O 2p bands. The Bi 6p and Pb 6p states, which oc- 
cur in the conduction bands from the CBM to ~ 4 eV 
above, also mix strongly with the O 2p bands. There is 
also cross-gap hybridization evident involving the Zn sp 
states, although this mixing is weaker perhaps because 
the Zn states occur further above the CBM. The values 
of the Born effective charges are enhanced by this cross 
gap hybridization. The average values for Bi, Pb, Zn and 
Nb are 4.4, 3.6, 2.9 and 5.8, respectively. The calculated 
polarization for the supercell is 85 /iC/cm 2 . This num- 
ber was obtained using a Berry's phase calculation. The 
estimate of the polarization that is obtained by multiply- 
ing the Born charges and cation off-centerings in their O 
cages is 83 fiC/cm 2 in close agreement with the Berry's 
phase calculation. For comparison the polarization of 
BiFe0 3 from DFT calculations is 90 - 100 /xC/cm 2 ^ 

The above discussion in terms of cross gap hybridiza- 



tion due to interaction between O 2p states and extended, 
low lying metal orbitals making up the conduction bands 
is different from a discussion in terms of a polarizable 
atomic like s orbital. However, it should be emphasized 
that atomic character is not a uniquely defined quantity 
since it is basis set dependent (note that an expansion 
in terms of radial functions times spherical harmonics 
about a single site is a complete basis, and so expan- 
sion about the different sites on a lattice is necessarily 
overcomplete) . Furthermore, polarization of an s orbital 
implies a mixing in of some p character from unoccu- 
pied orbitals, and in these perovskites this is 6p charac- 
ter that arises from conduction bands. The distinction 
between the two views is that in the lone pair description 
the discussion is in terms of response to crystal potential 
as in a discussion of crystal field, while in the cross-gap 
hybridization it is structure dependent hybridization as 
in ligand field. This is analogous to the origin of crys- 
tal field splittings in transition metal oxide compounds, 
which arise primarily from hybridization with O 2p states 
and not the non-spherical crystal potential around the 
metal site. 

As mentioned, high piezoelectric coupling in ferroelec- 
tric perovskites is generally found near an MPB, most 
commonly between a pseudocubic phase and a tetragonal 
phase, as in PZT. We did calculations as a function of im- 
posed tetragonal strain for our supercell at a fixed volume 
equal to the equilibrium volume for the pseudocubic cell 
(a=3.99 A). Over the range considered (0.992 < c/a < 
1.076), the cation off-centerings stayed collinear and pri- 
marily along a [111] direction, and in particular, at least 
up to the maximum strain considered, did not switch to 
an [001] average direction. As an example, cation off- 
centerings for c/a = 1.0606 are shown in Fig. [4] The 
lowest energy is at c/a ~ 1.015. However, near this min- 
imum the cation displacements remain very close to a 
[111] direction. Instead it seems that the deviation of c/a 
from unity is due to the octahedral tilts and not strain 
coupling for the tetragonal ferroelectric state. Since the 
tilts are expected to be coupled to the A-site cation or- 
dering, this may very well be an artifact of the specific 
ordering used in the necessarily small supercell studied 
here (note that the A-site cations would most likely be 
disordered in this material). In any case, BiPbZnNbOg is 
most likely pseudocubic with a rhombohedral ferroelec- 
tric state or weakly tetragonal. As such it may be of 
interest to investigate the BiPbZnNbOe ~ PbTi0 3 solid 
solution to determine whether an MPB is present. We 
note that the tetragonality may be particularly sensitive 
to volume, as is the case in PbTi0 3 , meaning that the 
LDA volume error may be particularly important for this, 
and also that errors due to the use of a small ordered su- 
percell may be particularly important for this. 

BiPbZnNbOe has two A-site ions with stereochemi- 
cal activity. Bi shows large off-centerings while the Pb 
off-centerings are smaller, and furthermore Bi carries a 
higher charge. This raises the question of the role of Pb 
in the ferroelectricity. This question is also of interest 
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TABLE II: Comparison of properties of BiPbZnNbOe and 
BiSrZnNbOg cells as obtained from our LDA calculations. 
a denotes the pseudocubic lattice parameter, the z* are the 
average Born effective charges, 8 are the average cation off- 
centerings with respect to the O cages, and P is the polariza- 
tion. 





BiPbZnNbOe 


BiSrZnNbOs 


a (LDA) 


3.99 A 


3.97 A 


*(Bi) 


0.64 A 


0.62 A 


**(Bi) 


4.4 


4.4 


<5(Pb/Sr) 


0.39 A 


0.33 A 


z*(Pb/Sr) 


3.6 


2.6 


<5(Zn) 


0.27 A 


0.25 A 


2*(Zn) 


2.9 


2.8 


<5(Nb) 


0.27 A 


0.25 A 


z*(Nb) 


5.8 


5.8 


P(lll) 


85 nC/cm 2 


79 /iC/cm 2 



because Pb is considered hazardous. This fact motivates 
searches for Pb-free replacements for the commonly used 
Pb-based ferroelectrics such as PZT. We performed cal- 
culations for 40 atom BiSrZnNbOg supercells, similar to 
those discussed above for the Pb system. Some properties 
are compared with the Pb system in Table [TT1 We found 
a slightly smaller equilibrium volume for the Sr system, 
corresponding to a pseudocubic a=3.97 A. As shown in 
Fig. [5j the material is ferroelectric, with large collinear 
cation off-centerings along [111]. The calculated polar- 
ization is 79 ^C/cm 2 . These off-centerings are similar to 
those for the Pb system, except on the Sr site. This is 
also the case for calculations at the same volume. For a 
pseudocubic a=3.98 A, the calculated off-centerings are 
0.62 A, 0.33 A, 0.25 A, and 0.25 A, for Bi, Sr, Zn and 
Nb, respectively, as compared with 0.62 A, 0.38 A, 0.26 
A, and 0.27 A, for Bi, Pb, Zn and Nb in the Pb compound 
at the same volume. This implies only slightly inferior 
ferroelectric properties for the Sr compound as compared 
to the Pb compound. The calculated average Born effec- 
tive charges for BiSrZnNb0 6 are 4.4, 2.6, 2.8, and 5.8 for 
Bi, Sr, Zn, and Nb, respectively. These are very similar 
to BiPbZnNbOg, except that the Sr Born charge is lower 
than that for Pb. Importantly, the LDA band gap of the 
Sr compound, £^=3.0 eV, is slightly higher than that of 
the Pb compound. This large band gap combined with 
the stable valences of the various cations implies that 



highly insulating samples can be prepared. This may 
make BiSrZnNbOg a useful alternative Pb-free ferroelec- 
tric. 



IV. SUMMARY AND DISCUSSION 

Density functional calculations for small supercells of 
double perovskite BiPbZnNbOg and BiSrZnNbOg indi- 
cate that these may be useful ferroelectric materials. 
They have high polarization due to the large displace- 
ment of Bi and the displacements of Nb combined with 
its high charge. They have substantial band gaps and 
stable chemistry consistent with synthesis of highly insu- 
lating samples. It will be of interest to measure prop- 
erties experimentally, and to investigate the solid so- 
lutions with PbTiOa to determine if a morphotropic 
phase boundary is present. Finally, we note that be- 
cause these materials have mixtures of ions with differ- 
ent charges on both the A- and -B-sites, they should 
have considerable chemical flexibility. While we con- 
sidered BiPbZnNbOg in the double perovskite structure, 
Bi x Pbi_ a: Zni/ 3+x /3Nb2/3_ a; /303 may exist over a large 
range of x. BiPbZnNbOg is the £=0.5, while x=Q is the 
known relaxor ferroelectric PZN. To our knowledge, the 
x=l perovskite BiZn 4/ / 3 Nbi/ 3 03 has not been reported, 
and instead the Bi203 - ZnO - Nb20s system shows Nb 
rich pyrochlore phases that presumably compete with the 
perovskite4i As such, it may be possible to synthesize the 
perovskite under pressure, similar to other low tolerance 
factor Bi based perovskites . 42 ' 43 In any case, the chemi- 
cal flexibility implies a tunability of properties that may 
be very helpful in obtaining useful ferroelectric properties 
including MPBs. 
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